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SUMMARY 

A  unified  set  of  composite  micromechanics  equations  of  simple  form  Is 
.summarized  and  described.  This  unified  set  Includes  composite  micromechanics 
equations  for  predicting  (1)  ply  In-plane  uniaxial  strengths;  (2)  through-the- 
thlckness  strength  (Interlaminar  and  flexural);  (3)  In-plane  fracture  tough¬ 
ness;  (4)  In-plane  Impact  resistance;  and  (5)  through-the-thlckness  (Inter¬ 
laminar  and  flexural)  Impact  resistance.  Equations  are  also  Included  for 
predicting  the  hygrothermal  effects  on  strength,  fracture  toughness  and  Impact 
resistance.  Several  numerical  examples  are  worked  out  to  Illustrate  the  ease 
of  use  of  the  various  composite  micromechanics  equations.  The  numerical 
examples  were  selected,  In  part,  to  demonstrate  the  Interrelationships  of  the 
various  constituent  properties  In  composite  strength  and  strength-related 
behavior,  to  make  comparisons  with  available  experimental  data  and  to  provide 
Insight  Into  composite  strength  behavior. 


INTRODUCTION 

The  several  strengths  (stresses  at  fracture)  of  unidirectional  composites 
are  fundamental  to  analysis/design  of  fiber  composite  structures.  Some  of 
these  strengths  are  determined  by  physical  experiments.  Others  are  not  easily 
amenable  to  direct  measurement  by  testing.  In  addition,  testing  Is  usually 
time  consuming,  costly,  and  the  composite  must  have  been  made  prior  to  testing. 
Furthermore,  parametric  studies  of  the  effects  of  fiber  volume  ratio  on  proper¬ 
ties  such  as  Impact  resistance  and  fracture  tojughness  can  only  be  made  by  an 
extensive  combination  of  tests.  Another  appxbach  Is  the  use  of  composite 
micromechanics  to  derive  equations  for  predicting  composite  strengths  based  on 
constituent  (fiber  and  matrix)  properties.  Over  the  last  twenty  years,  com¬ 
posite  micromechanics  has  been  used  to  derive  equations  for  predicting  selected 
composite  strengths  (ref.-.l).  However,  these  equations  are  not  readily  avail¬ 
able  since  equations  for  different  strengths  are  scattered  throughout  the 
literature. 

Herein,  a  unified  set  of  composite  micromechanics  equations  Is  summarized 
and  described.  The  set  Includes  simple  equations  for  predicting  ply  (uni¬ 
directional  composite)  strengths  using  constituent  properties.  Equations  are 
for:  (1)  tensile  strengths  (In-plane  and  through-the-thlckness),  (2)  flexural 
strength,  (3)  Impact  resistance,  and  (4.)  fracture  toughness.  Also,  equations 
are  presented  for  the  effects  of  (1)  moisture,  and  (2)  temperature.  Results 
predicted  by  these  equations  are  compared  with  available  experimental  data. 
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These  data  are  primarily  from  Refs.  1  to  5.  The  equations  are  summarized  In 
subsets  corresponding  to  related  strengths  such  as  In-plane,  through-the- 
thlckness,  etc.  The  description  consists  of  the  significance  of  the  partici¬ 
pating  variables  In  the  equations  of  each  subset,  several  numerical  examples 
and  possible  Implications. 

The  equations  of  each  subset  (strengths,  fracture  toughness,  Impact  resis¬ 
tance  and  hygorthermal  degradation  effects)  are  summarized  In  chart  form 
(labeled  figures).  This  allows  the  equations  for  each  subset  to  be  In  one  page 
for  convenience  of  use  and  Identification  of  Interrelationships.  Constituent 
material  properties  used  In  the  numerical  examples  are  tabulated  and  Identified 
with  the  same  symbol  used  In  the  equations.  The  numerical  examples  are  pre¬ 
sented  In  narrative  form,  rather  than  tabular,  In  order  to  conserve  space.  The 
symbols  used  are  summarized  In  the  Appendix  for  convenience  of  reference. 

Many  of  the  equations  Included  In  this  composite  micromechanics  unified 
set  appear  In  their  present  simplified  form  for  the  first  time.  These  equa¬ 
tions  evolved  from  continuing  research  on  composite  micromechanics  and  com¬ 
posite  computational  mechanics  at  Lewis  Research  Center.  Also,  this  Is  the 
first  unified  set  which  provides  a  quantified  description  of  composite  strength 
and  strength-related  behavior  (fracture  toughness  and  Impact  resistance)  at  the 
mlcromechanlstlc  level. 


SYMBOLS 

c  heat  capacity 

D  dlffuslvlty 

d  diameter 

E  modulus  of  elasticity 

G  shear  modulus 

Impact  energy  density 
K  heat  conductivity 

k  volume  ratio 

M  moisture  -  percent  by  weight 

Nf  number  of  filaments  per  roving  end 

P  property 

RHR  relative  humidity  ratio 

S  strength 

d.  fracture  "toughness" 

T  temperature 

t  thickness 

x,y,z  structural  reference  axes 

1,2,3  ply  material  axes 

a  thermal  expansion  coefficient 

B  moisture  expansion  coefficient 
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<5  Interfiber,  Interply  spacing 

5  fracture  strain,  strain 

e  ply  orientation  angle 

X  weight  percent 

p  density 

a  stress 

Subscripts: 

F  fiber  property 

C  compression  property 

0  dry  property 

6  glass-transition 

F  flexural 

8.  ply  property 

m  matrix  property 

S  shear 

SB  short  beam  shear 

T  tension 

v  void 

W  wet 

o  reference  property,  temperature 

«  saturation 

1,2,3  direction  corresponding  to  1,2,3  ply  material  axes 


COMPOSITE  MECHANICS  —  DEFINITIONS  AND  CONSTITUENT  MATERIALS 

The  branch  of  composite  mechanics  which  provides  the  formal  structure  to 
relate  ply  uniaxial  strengths  to  constituent  properties  Is  called  composite 
micromechanics.  Composite  micromechanics  for  uniaxial  strengths  Is  Identified 
concisely  In  the  schematic  In  figure  1.  The  schematic  In  this  figure  defines 
the  Inputs  to  composite  micromechanics  and  the  outputs.  The  Inputs  consist  of 
constituent  material  (fiber/matrix)  properties,  geometric  configuration,  en¬ 
vironmental  condition's,  and  the  fabrication  process.  The  outputs  consist  of 
ply  uniaxial  strengths,  Impact  resistance,  fracture  toughness  and  hygrothermal 
effects. 

The  formal  structure  of  composite  micromechanics  (concepts,  math-models 
and  equations)  Is  developed  based  on  certain  assumptions  (consistent  with  the 
physical  situation)  and  the  principles  of-  solid  mechanics.  The  four  main 
assumptions  made  In  deriving  the  equations  described  herein  are:  (1)  the  ply 
resists  loads  as  depicted  schematically  In  figure  2;  (2)  the  ply  and  Its  con¬ 
stituents  behave  linearly  elastic  to  fracture  as  Is  Illustrated  In  figure  3; 
(3)  the  ply  uniaxial  strengths  are  associated  with  their  respective  fracture 
modes  shown  In  figure  4;  and  (4)  there  Is  complete  bond  at  the  Interface  of 
the  constituents.  Though  the  principles  of  solid  mechanics  can  be  used  with 
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various  levels  of  mathematical  sophistication,  the  mechanics  of  materials  was 
used  In  deriving  the  equations  summarized  herein  because  It  leads  to  explicit 
equations  of  simple  form  for  each  property. 

Properties  along  the  fiber  direction  (1-axls,  fig.  2)  are  conventionally 
called  longitudinal;  those  transverse  to  the  fiber  direction  (2-axls,  fig.  2) 
are  called  transverse;  the  In-plane  shear  Is  also  called  Intralaminar  shear 
(1-2  plane,  fig.  2).  Those  through  the  thickness  (3-axls,  fig.  2)  are  called 
Interlaminar  properties.  All  ply  properties  are  defined  with  respect  to  the 
ply  material  axes  denoted  by  1 ,  2  and  3  In  figure  2  for  descrlptlon/analysls 
purposes.  Most  ply  properties  are  denoted  by  a  letter  with  suitable  sub¬ 
scripts.  The  subscripts  are  selected  to  Identify  type  of  property  (ply,  fiber, 
matrix),  plane,  direction,  and  sense  In  the  case  of  strengths.  For  example, 
SfclU  denotes  ply  longitudinal  tensile  strength  while  Sfj  denotes  fiber  tensile 
strength.  Though  this  notation  may  seem  cumbersome.  It  Is  necessary  to  prop¬ 
erly  differentiate  among  the  multitude  of  ply  and  constituent  properties. 

A  variety  of  fibers  have  been  used  to  make  composites.  Some  of  these  are 
summarized  In  table  1  with  their  respective  properties  needed  for  composite 
micromechanics.  Similarly,  some  typical  matrix  resins  are  summarized  In 
table  2. 


UNIAXIAL  STRENGTHS  —  IN-PLANE 

There  are  five  In-plane  ply  uniaxial  strengths.  These  are  Identified  as: 
(1)  longitudinal  tension  (S^nj);  (2)  longitudinal  compression  (Sj^q);  (3) 
transverse  tension  (S^t);  (4)  transverse  compression  (Sj^c);  a°d  (5)  In¬ 
plane  or  Intralaminar  shear  (Sa-|2s)*  The  fracture  modes  associated  with  each 
uniaxial  strength  are  depicted  schematically  In  figure  4.  Note  that  there  are 
three  different  and  distinct  fracture  modes  for  longitudinal  compression 
(fig.  4-b) :  (1)  fiber  compression  (shear  plane)  fracture;  (2)  delamlnatlon 

transverse  splitting  or  panel  buckling;  and  (3)  fiber  microbuck  ling. 

The  composite  micromechanics  equations  for  the  ply  uniaxial  strengths  are 
summarized  In  figure  5  with  attendant  schematics.  The  schematics  define  the 
load  direction,  fiber  orientation  and  the  notation  used  In  the  micromechanics 
equations.  The  first  five  equations  describe  the  In-plane  uniaxial  strengths, 
respectively:  S^nf,  Sj,nc>  %22T»  S  8.22C*  and  $sn2S*  The  ^as^  equation  Is 
for  the  void  effect  on  the  resin  strength  (Sm)  and  also  provides  lower  bound 
estimates  on  S^22T»  %22C  ar,d  S112S  as  1)6  described  later. 

The  following  are  observed  from  the  micromechanics  equations  for  ply  uni¬ 
axial  strengths:  (T)  depends  on  Sfj  and  th6  fiber  compression  fracture 

mode  for  Sanc  depends  on  SfC.  These  are  the  only  two  that  are  fiber  strength 
dominated.  (2)  The  delamlnatlon/spllttlng  for  S^c  depends  on  matrix  shear 
strength  (through  Sa-|2S  (e(l*  5))  and  the  matrix  tensile  strength  and,  there¬ 
fore,  Is  resin  strength  dominated.  (3)  The  mlcrobuckllng  fracture  mode  for 
SilIC  depends  strongly  on  the  shear  modulus  and  mildly  on  the  modular  ratio 
(Gm/Gf|2)  and,  therefore,  Is  resin  stiffness  dominated.  (4)  The  other  three 
(si22T»  $122C  a°d  sn2s)  depend  strongly  on  the  respective  resin  strengths  and 
are,  therefore,  resin  strength  dominated.  (5)  The  fiber- volume  ratio  affects 
strongly  and  Sj^-jc  (fiber  compressing  or  fiber  mlcrobuckllng)  which  are 

the  fiber  strength  or  shear  stiffness  dominated  ply  strengths).  (6)  The  fiber 
volume  ratio  affects  mildly  Sj^t*  sl22T*  S!22C»  S112S  and  S111C  (delamlnatlon/ 
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shear)  which  are  the  resin-strength  dominated  ply  uniaxial  strengths.  (7)  The 
voids  Influence  the  matrix  strength.  Several  examples  below  Illustrate  use  of 
the  uniaxial  strength  equations  In  figure  5. 

Example  3.1.  Calculate  the  ply  tensile  strength  (Sty|-|-)  of  a  graphite 
fiber  (AS)/1ntermed1ate  modulus,  high  strength  (IMHS)  epoxy  (AS/E)  com¬ 
posite  with  0.6  fiber  volume  ratio.  From  table  1,  Sft  =  350  ks 1  and  from 
equation  (1)  (fig.  5),  S^ny  =  210  ks  1  which  Is  the  same  as  the  measured 
value  In  table  3. 

Example  3.2.  Calculate  the  ply  compression  strength  for  an  AS/E  (IMHS) 
composite  with  0.6  fiber  volume  ratio.  All  three  equations  (eq.  (2), 
fig.  5)  should  be  used.  To  use  the  first  equation,  Sfc  must  be  known.  If 
It  Is  not  known-,  SfC  «  0.9  Syy  Is  a  good  approximation  for  graphite  fiber/ 
matrix  composites  (ref.  6).  Using  this  approximation  and  respective 
values  for  tables  1  and  2  In  equation  (2)  figure  5  from  the  first 
equation: 

SftHC  -  0.6  x  °-9  x  350  =  189  ks  1 

from  the  second  equation  (note  need  evaluate  Sjh2S  from  eq.  (5),  example 
3.5  with  Incomplete  bond) 

S411C  =  10.. 0  x  8.1  +  2.5  x  15  =  118  ksl 
from  the  third  equation, 

SjnC  *  0.185/[1 .0-0.6  (1.0-0.185/2.0)]  =  406  ksl 

A  conservative  approach  Is  to  select  the  lowest  value  or  S^hq  =  118  ksl. 
This  Is  about  69  percent  of  the  typical  measured  value  of  170  ksl  In  table  3. 
The  value  of  189  ksl  predicted  by  the  first  equation  Is  also  reasonable.  This 
value  could  be  used  In  laminates  which  have  other  than  0°  plies  on  the  outside. 

Example  3.3.  Calculate  the  ply  transverse  tensile  strength  (Sj^t)  for 
an  AS/E  (IMHS)  composite  with  0.6  fiber  volume  ratio.  From  equation  (3) 
and  respective  values  from  tables  1  and  2 
S422T  »  [1.0  -  ( -y/0.6  -  0.6)  x  (1.0  -  0. 5/2.0)]  x  15  =  13  ksl 

This  Is  about  twice  the  measured  value  of  7  ksl  In  table  3.  One  major 
reason  for  this  difference  Is  the  complete  bond  at  the  fiber/matrix  Interface 
assumed  In  deriving  equation  (3),  figure  5.  .Incomplete  bond  at  the  Interface 
may  be  approximated  by  assuming  the  presence  of  voids.  Assuming  about  5 
percent  voids  by  volume  (kv  =  0.05)  and  using  equation  (6)  (fig.  5),  the 
reduced  or  degraded  resin  tensile  strength 

Smy  -  {l.O  -  [4  (0.05) /(1 .0-0.6)  ir]l/?}  x  15  =  9.0  ksl 

Using  this  reduced  resin  strength  In  equation'  (3)  (fig.  5)  S^22T  *  7-8  ksl 
which  Is  a  reasonable  estimation  compared  to  the  measured  value  of  7  ksl. 

The  above  calculations  lead  to  the  conclusion  that  estimation  of  Sa22T> 
then,  requires  two  steps: 

1.  Degradation  of  Smy  due  to  5  percent  voids  by  volume  (kv  =  0.05)  as 
predicted  by  equation  6  (fig.  5). 

2.  Substitution  of  the  degraded  Smy  In  equation  (3j. 

Though  the  value  of  5  percent  voids  may  seem  somewhat  arbitrary,  It  Is 
reasonable  since  equation  (3),  figure  5,  does  not  account  for  factors  such  as 
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nonuniform  fiber  distribution  within  the  ply,  Incomplete  (partial)  bond  at  the 
fiber/matrix  Interface,  and  possible  differences  In  the  In  situ  resin  matrix 
properties  compared  to  the  neat  resin  properties. 

Example  3.4.  Calculate  the  transverse  compression  strength  (Sa22C)  for 
an  AS/E  ( IMHS)  with  0.6  fiber  volume  ratio.  From  equation  (4)  and  respec¬ 
tive  values  from  tables  1  and  2 

S122C  =  H-0  "  (^0  -  0.6)  (1.0  -  0. 5/2.0)]  x  35.0  =  30.4  ksl 

This  values  Is  about  85  percent  of  the  typical  value  In  table  3.  It  Is 
worth  noting  that  the  Interfaclal  bond  and  fiber  nonuniformity  are  not  critical 
In  transverse  compression  and,  therefore,  do  not  contribute  to  resin  compres¬ 
sion  strength  degradation. 

Example  3.5.  Calculate  the  Intralaminar  shear  strength  (Sj^s)  f°r  an 
AS/E  (IMHS)  with  0.6  fiber  volume  ratio.  From  equation  (5)  and  respective 
values  from  tables  1  and  2 

Sj_i2S  -  [1.0  -  (/O  -  0.6)  (1.0  -  0.185/2.0)  ]xl3  =  10.9  ksl 

This  values  Is  about  21  percent  greater  than  the  typical  measured  value 
9  ksl,  table  3.  This  value  Is  reasonable  In  view  of  the  nonuniform  fiber  dis¬ 
tribution  bonding  condition  at  the  Interface  and  deviations  In  In  situ  proper¬ 
ties  from  neat  resin  properties  as  was  mentioned  for  S&22T  (example  3.2).  A 
closer  estimate  to  the  measured  value  may  be  obtained  by  degrading  the  matrix 
shear  strength  Sms  assuming  2  percent  voids  (kv  =  0.02)  In  equation  (6), 
figure  5.  The  result  Is  Sj^s  =  8.1.  It  Is  worth  noting  that  the  various  fac¬ 
tors  that  affect  the  transverse  tensile  strength  also  affect  the  Intralaminar 
shear  strength  but  not  as  severely. 

Example  3.6.  Calculate  the  effect  of  voids  on  the  ply  transverse  tensile 
strength  of  an  AS/E  (IMHS)  with  0.6  fiber  volume  ratio  (kf  =  0.6)  and  0.02 
void  volume  ratio  (kv  =  0.02).  This  Is  accomplished  using  the  following 
three  steps: 

1.  Voids  effect  on  Smj  (eq.  (6),  fig.  5). 

SmT  =  {l.O  -  [4(0.02)/(l .0  -  0.6)  x  15  =  11.2  ksl 

2.  Incomplete  Interfaclal  bond  effects' 

SmT  =  |1.0  -  [4  (0.05)/(l  .0  -  0.6  ir]1/2i  x  11.2  =  6.8  ksl 

3.  Transverse  tensile  strength  (eq.  (3),  fig.  5) 

S122T  =  [1.0  -  (yO“-  0.6)  (1  -  0. 5/2.0)]  6.8  =  5.9  ksl 

Note  that  the  two  void  ratios  (kv  =  0.02  and  kv  =  0.05)  are  not  additive 
since  equation  (6)  Is  nonlinear. 

Example  3.7.  Calculate  the  lower  bound  of  the  ply  transverse  tensile 
strength  for  an  AS/E  (IMHS)  composite  with  kf  =  0.6.  For  this  case,  we 
use  equation  (6),  figure  5,  with:  kv  =  0.6,  kf  =  0.0  and  S^T  =  %• 

S&22T  =  {1.0  -  [4  (0.6)/ir]l/2|  x  15  =  1.9  ksl 

This  value  Is. about  27  percent  of  the  typical  measured  value  of  7  "ksl  and 
about  15  percent  of  that  predicted  by  equation  (3),  figure  5,  without  any  deg¬ 
radation.  The  lower  bound  Is  overly  pessimistic  for  acceptable  composites  and 
should  be  used  only  by  In .composites  with  no  Interfaclal  bond.  Lower  bound 
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estimates  on  ply  transverse  compression  (Sj?2c)  and  ply  Intralaminar  shear 
strength  (Sj^s)  are  obtained  by  following  the  same  procedure. 


UNIAXIAL  STRENGTHS  —  THROUGH-THE-THICKNESS 

There  are  six  through-the-thlckness  uniaxial  strengths.  These  are  Iden¬ 
tified  as:  (1)  longitudinal  Interlaminar  shear  (parallel  to  the  fiber  direc¬ 
tion),  (Sji3$);  (2)  transverse  Interlaminar  shear  (transverse  to  the  fiber 
direction),  (Sa23s);  (3)  longitudinal  short-beam-shear  (parallel  to  the  fiber 
direction),  (Sj^sb);  (4)  transverse  short-beam-shear  (transverse  to  the  fiber 
direction),  (S123sb);  (5)  longitudinal  flexural  (bending)  (Sayjp);  and  (6) 
transverse  flexural  (S^22p).  The  composite  micromechanics  equations  for  these 
uniaxial  strengths  are  summarized  In  figure  6  with  attendant  schematics.  The 
first  six  equations  describe  the  six  through-the-thlckness  uniaxial  strengths, 

respectively:  S113c,  Sa23S.  %13SB-  st23SB»  slllF  and  %22F-  The  ]astu 
equation  describes  the  void  effects  on  the  resin  strength  and  can  also  be  used 

as  a  lower  bound  on  ply  strengths  dominated  by  the  resin  as  was  mentioned 
previously. 

The  following  are  observed  from  the  composite  micromechanics  equations  In 
figure  6:  (1)  Sj^os  Is  the  same  as  Sj^2s;  (2)  S423S  depends  strongly  on  the 

resin  shear  strength  (Sm$)  and  mildly  on  kf  and  Gm/Gf23i  (3)  the  short-beam- 
shear  strengths  S^sb  and  S^jb  are  1.5  times  their  respective  Interlaminar 
shear  strengths  ( Sa-] 05  and  Sa23s);  (4)  the  longitudinal  flexural  strength 
( Sft,i ip)  fiber  dominated  ana,  thus,  depends  strongly  on  kf,  spy  and  Sf^;  (5) 
the  transverse  flexural  strength  (Sa22p)  Is  matrix  strength  dominated  and,  thus, 
depends  strongly  on  SmT  and  Smc  but  It  depends  mildly  on  kf  and  Em/Ef22;  (6) 
the  voids  degrade  matrix  strength  depending  nonllnearly  on  both  kv  ana  kf. 
Several  examples  below  Illustrate  use  of  the  equations  In  figure  6. 


Example  4.1.  Calculate  the  longitudinal  Interlaminar  shear  strength 
(SjH3$),  for  an  AS/E  (IMHS)  composite  with  0.6  fiber  volume  ratio.  This 
Is  the  same  as  Sa-|2s.  However,  we  go  through  the  steps  again  for  com¬ 
pleteness.  Using  equation  (4),  the  first  of  equation  (1),  figure  6,  and 
respective  property  values  from  tables  1  -and  2: 


1.  Incomplete  bond  simulation  (kv  =  0.02) 

smS  -  l1-0  "  t4  ( 0 . 02 ) /( 1 .0  -  0.6)ir]1/2}  x  13  =  9.7  ksl 

2.  Longitudinal  Interlaminar  shear  strength 

=  n*°  -  -  0.6)  0  -  0.185/2.0)]  x  9.7  =  8.1  ksi 

which  Is  the  sameT as  Sa-|2s  as  expected. 


Example  4.2.  Calculate  the  transverse  Interlaminar  shear  strength  (St235) 
for  an  AS/E  (IMHS)  with  0.6  fiber  volume  ratio.  Using  Sm5  from  example 
4.1,  figure  6,  and  respective  properties  from  table  1  and  2  In  the  second 
equation  (1),  figure  6: 


p.O  -  t/Q~  (1.0  -  0.185/1.0  1 
H23S  =  Ll.O  -  0.6  (1.0  -  0. 185/1. 0.)J  X 


7.0  ksl 


Note  that  Sj_23s  0.85  Say 35  Indicating  that  the  ply  Is  weaker  In 
transverse  Interlaminar  shear  than  In  longitudinal  Interlaminar  shear  strength. 
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Example  4.3.  Calculate  (1)  the  longitudinal  short-beam-shear  ( 3S8)  f°r 
the  composite  In  example  4.1;  and  (2)  the  transverse  short-beam-shear 
(S823SB>  1n  example  4.2. 

1.  Using  Sfci3s  _  g.i  in  the  first  of  equation  (2),  figure  6 
s!ll 3SB  =  ]-5  x  8.1  =  12.1  ks  1 

This  value  Is  In  reasonably  good  agreement,  under  estimating  by  14 
percent,  with  the  typical  measured  value  of  14  ks 1  In  table  3. 

2.  Using  Sj^23S  -  7.0  In  the  second  of  equation  (2),  figure  6 
%23SB  =  1*5  x  7.0  =  10.5  ksl 

Measured  values  for  this  strength  are  not  available  for  comparison. 


Example  4.4.  Calculate  the  ply  longitudinal  flexural  strength  (Soyir)  for 
an  AS/E  (IMHS)  composite  with  0.6  fiber  volume  ratio.  Using  kf  =  0.6, 

Sfc  =  0.9  Sfj,  and  Sf-j-  *  350  ksl  (table  1)  in  the  first  of  equation  (3), 
figure  6. 


.  3  (0.61  [ 350) 

*8.1 1 F  “  1.0  +  1.0/0. 9 


=  298  ksl 


This  value  over  estimates  the  typical  measured  data  of  230  ksl  by  about 
30  percent.  A  lower  estimate  for  this  strength  Is  obtained  by  using 
SallT  =  210  ksl  (from  example  3.1)  and  Sollc  =  118  ksl  (from  example  3.2)  In 
the  following  equation: 


3  S 


81 1 T 


81  IF 


1  +  S81 1 T/S8l 1 C 


3  X  210 
1.0  +  210/118 


227  ksl 


which  estimates  the  measure  value  of  230  ksl  almost  exactly.  This  Is  defi¬ 
nitely  a  very  good  estimate  considering  the  simplicity  of  the  equations  and  the 
uncertainties  associated  with  longitudinal  compression  failure  (ref.  6).  It 
also  Illustrates,  In  part,  that  flexural  failure  Is  probably  a  complex  combi¬ 
nation  of  tension,  compression,  and  Intralaminar  shear  failures. 


15  =  27.4  ksl 


Example  4.5.  Calculate  the  ply  transverse  flexural  strength  of  an  AS/E 
(IMHS)  composite  with  0.6  fiber  volume  ratio.  Using  kf  =  0.6  and  re¬ 
spective  property  values  from  tables  1  and  2  In  the  second  of  equation  3. 
c  j 3  [1.0  -  (t/O"-  0.6)  (1.0  -  0.5/2.01  * 

822F  j 1.0  +15/35 

This  over  estimates  the  typical  value  of  18  ksl  (table  3)  by  about  52  percent. 
As  was  the  case  for  longitudinal  flexural  strength,  a  lower  estimate  can  be 
obtained  by  substituting  Sjt22T  =  7-8  ksl  (from  example  3.3  with  partial  Inter¬ 
facial  bond)  and  S  a,22C  ~  80.4  ksl  (example  3.4)  In,  the  following  .equation 

S  3  X  S822T  3  x  7.8  .  . 

822F  1  +•  Sa.22T/Sjt22C  "  1,0  +  7-8/30*4  =  18'6  ks1 


which  Is  almost  equal  to  the  typical  measured  values  of  18  ksl  (table  3). 

Example  4.6.  Calculate  the  effect  of’  3  percent  voids  on  the  longitudinal 
flexural  strength  for  an  AS/E  (IMHS)  with  0.6  fiber  volume  ratio.  The 
first  of  equation  (3)  shows  no  void  effect.  However,  the  lower  estimate 
equation  In  example  (4.4)  Indicated  that  the  compression  strength  pre¬ 
dicted  by  the  second  of  equation  (3)  should  be  used.  This  Is  calculated 
using  the  following  steps: 


8 


1. 

2. 


3. 


4. 

5. 


Sms  degraded 
SmS  -  {1.0  - 
Intralaminar 
step  1 . 

S112S  =  P  - 
Smj  degraded 

smS  , 

SmT  =  {1.0  - 
Longitudinal 
Sine  =  TO-0 
Longitudinal 


for  voids  (eq.  (4),  fig.  6) 

[4  ( 0 . 03 ) /( 1 .0  -  0.6)  w]1/*)  x  13  =  9.0  ksl 

shear  strength  Sa-|2$  (eq.  (5),  fig.  5)  with  Sm$  from 

(  0.6  -  0.6)  (1.0  -  0.185/2.0)]  9.0  =  7.6  ksl 

for  voids,  void  degradation  ratio  same  as  In  step  1  for 

[4  (0.03)/(l  .0  -  0.6)  ir]!/2}  x  15  =  10.4  ksl 
compression  S^ne  (second  of  eq.  (2),  fig.  5) 
x  7.6  +  1.5  x  10.4  =  102  ksl 

flexural  strength  (lower  estimate  equation,  example  4.4) 


SlllF  =  3  x  210/(1.0  +  210/102)  =  206  ksl 


which  Is  about  9  percent  less  than  the  227  ksl  value  calculated  without 
voids  In  example  4.4.  Two  points  are  worth  noting:  (1)  step  (4)  results 
In  a  void  degradation  of  about  13  percent  In  Sayic  and  (2)  the  void  deg- 
radatlon  Is  more  severe  for  the  longitudinal  compression  strength  than 
for  the  longitudinal  flexural  strength  (13  percent  versus  9  percent 
respectively) . 


The  above  calculations  show  that  the  composite  micromechanics  equations 
In  figure  6  and  the  alternates  In  Examples  4.4  and  4.5  can  be  used  to  obtain 
reasonable  estimates  for  through-the-thlckness  uniaxial  strengths.  The  calcu¬ 
lations  also  show  that  the  equations  can  be  used  to  Interpret  measured  data. 

In  either  case,  these  equations  should  be  used  judiciously. 


UNIAXIAL  FRACTURE  "TOUGHNESS" 

Fracture  toughness  Is  a  measure  of  a  material  to  resist  defects  such  as 
holes,  slits  and  notches.  Fracture  toughness  Is  described  by  fracture  tough¬ 
ness  parameters  associated  with  distinct  fracture  modes  (ref.  7).  Three  frac¬ 
ture  modes  are  generally  considered:  opening  mode  (Mode  I),  In-plane  shear 
(Mode  II)  and  out-of-plane  shear  (Mode  III).  In  the  case  of  unidirectional 
composites  and  assuming  full  thickness  penetration  defects,  there  are  three 
major  In-plane  fracture  toughness  parameters  herein  defined  as:  (1)  longitud¬ 
inal  fracture  toughness  Ce/niT);  (2)  transve.fse  fracture  toughness  le/j.22 t); 
and  (3)  In-plane  (Intralaminar)  shear  fracture  toughness  Ce/ii2S)-  These  three 
are  parallel  to  the  uniaxial  In-plane  strengths  S^yiy,  Si22f  and  sll2S*  respec¬ 
tively.  These  fracture  toughness  parameters  are  used  herein  to  denote  far- 
f 1  eld  stress  required  to  -produce  additional  damage  In  the  composite.  It  Is  not 
clear  whether  far  field  shear  stress  will  produce  Mode  II  fracture  In  unidirec¬ 
tional  composites  or  some  component  of  Mode  I  (opening)  'fracture.  In  view  of 
this  we  consider  only  Mode  I,  opening  fracture  modes  my  and  122T’ 

The  equations  describing  the  longitudinal  and  transverse  fracture  tough¬ 
ness  parameters  CefjniT  and  -®122T)  are  given  In  figure  7  (ref.  8)  with  attend¬ 
ant  schematics.  Two  sets  of  equations  are  given.  In  the  first  set  (eq.  (1) 
and  (2)),  mT  and  122T  are  expressed  In  terms  of  ply  properties  while  In 
the  second  set  (eq.  (3)  and  (4))  they  are  expressed  In  terms  of  constituent  properties. 
It  can  be  seen  In  Eqs.  (3)  and  (4)  that:  (1)  -e/mj  depends  linearly  on  Sf y  and 
depends  In  a  complex  way  on  Kf,  (Efll/Em),  (Em/Ef22),  <Gm/Gfl2)  and  vm;  and 
(2)  jd  122T  depends  linearly  on  SmT  and  In  a  complex  way  on  the  other  constit¬ 
uent  material  properties.  Examining  only  key  parameters  kf,  Ef  yj ,  Em,  Sf j  and 
SmT*  ^etlllT  Increases  with  Increasing  kf,  Efn ,  Em  and  with  decreasing  Ef  yj , 
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while  Increases  with  Increasing  Efn,  SmT  and  with  decreasing  kf  and  Em. 

Equations  (3)  and  (4),  figure  7,  are  cumbersome  to  use.  It  Is  easier  to  use 
equation  (1)  and  (2),  figure  7,  In  conjunction  with  the  ply  mechanical  proper¬ 
ties  equations  summarized  In  figure  8.  Two  examples  below  Illustrate  use  of 
these  equations  and  Interpretation  of  the  results. 

Example  5.1.  Calculate  the  longitudinal  fracture  toughness  of  an  AS/E 
(IMHS)  composite  with  0.6  fiber  volume  ratio.  This  Is  accomplished  using 
the  following  steps  together  with  equations  from  figures  7  and  8,  and 
respective  constituent  properties  from  tables  1  and  2. 


1.  Calculate  S^y  (eQ-  0).  fig-  5) 

SallT  =  0.6  x  350  =  210  ks 1 

2.  Calculate  Ejtn  (first  equation,  fig.  8) 

Em  =  0.6  x  31  +  0.4  x  0.5  =  18.8  mpsl 

3.  Calculate  E^22  (second  equation,  fig.  8) 

Ej?2  =  0. 5/[l  .0  -yO(l.  0-0. 5/2.0)]  =  1.2  mpsl 

4.  Calculate  »-)2  ( f Ifxh  equation,  fig.  8) 
vll 2  =  0*6  x  0.20  +  0.4  x  0.35  =  0.26 

5.  Calculate  G^i 2  (third  equation,  fig.  8) 

Gj^-j g  =  0.185/[1. 0  --i/076(  1-0. 185/2.0)]  =  0.62  mpsl 

6.  Calculate  jzfmj  (eq-  O).  fig-  7) 

^/oilT  =  210/  fl.O  +  [2  (18.8/1.2  -  0.26)  +  ( 18.8/0. 62)  j1 /2{ 

=  24  ks 1 


This  Implies  that  the  far-fleld  (P/A-type)  ply  stress  will  be  24  ks 1  when 
the  crack-like  defect  starts  growing.  This  can  also  be  Interpreted  as  follows: 
The  stress  required  to  produce  additional  damage  Is  reduced  by  a  factor  of 
about  ten  compared  to  that  In  a  ply  without  defects. 


Example  5.2.  Calculate  the  transverse  fracture  toughness  (-e/j,22T)  1°  an 
AS/E  (IMHS)  composite  with  0.6  fiber  volume  ratio.  This  Is  calculated 
using  appropriate  equations  figures  7  and  8,  respective  properties  from 
Tables  1  and  2  and  the  following  steps. 

1.  Calculate  Egj-|.  From  example  5.1  step  2,  E^n  =  18-8  mPs1 

2.  Calculate  Sa22T  (eq.  (3),  fl9-  5)  w1*h  partial  Interfacial  bond. 

a.  SmT  =  jl  .0  -  [4  (0.05)/(l  .0  ^0.6)  ir]1/2|  15  =  9.0ks1 

b.  S111T  =  [1-0  -  (/ITF-  0.6)  (f.O  -  0. 5/2.0)]  x  9.0  =  7.8  ksl 

3.  Calculate  Ejl22 *  Fro™  example  5.1,  step  3,  E^22  =  1-2  mpsl 

4.  Calculate  vj^2-  ^rom  example  5.1,  step  4,  vjn2  =  0.26 

5.  Calculate  Goi2*.  from  example  5.1,  step  5,  G^-j 2  =  0.62  mpsl 

6.  Calculate  >e3si22T-  From  equation  (2),  figure  7 


2.22T 


=  7.8/ 


0  +  (1.2/18. 8)1/2  [2  (1.0  -  0.26)' 


18.8/0.62]1/2}  =  3.2  ksl 


This  Implies  that  the  far-fleld  (P/A-type)  ply  stress  will  be  3.2  ksl  when 
additional  damage  In  the  vicinity  of  the  defect  will  occur.  Or  alternatively, 
the  stress  required  to  produce  additional  damage  Is  reduced  by  a  factor  of 
about  2.5  compared  to  that  In  a  ply  without  defects.  It  Is  worth  noting  that 
this  relatively  low  value  of  3.2  ksl  required  to  produce  additional  damage  Is 
a  major  contributor  to  the  brittle-like  strength  behavior-  transverse  to  the 
fiber  direction. 
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UNIAXIAL  IMPACT  RESISTANCE- IN-PLANE 

Uniaxial  Impact  resistance  of  unidirectional  composites  Is  defined  herein 
as  an  In-plane  uniaxial  Impact  energy  density.  It  Is  denoted  by  the  generic 
symbol  and  Is  associated  with  a  corresponding  In-plane  uniaxial  Impact 
stress.  There  are  five  Impact  energy  densities:  (1)  longitudinal  tension 
(jE^ht);  (2)  longitudinal  compression  (JE^nc);  (3)  transverse  tension  (j£&22T); 

(4)  transverse  compression  ( %%220  and  Intralaminar  shear  (#stl2s)-  The  composite 
micromechanics  equations  for  these  Impact  energy  densities  are  summarized  In 
figure  9  with  attendant  schematics.  The  wlggly  arrows  In  the  schematics  denote 
dynamic  stresses.  These  equations  are  derived  by  assuming  linear  stress-strain 
behavior  to  fracture  (fig.  3)  under  dynamic  stress.  The  first  five  equations 
describe  the  five  In-plane  uniaxial  Impact  energy  densities  while  the  last 
equation  describe  the  void  degradation  effect  as  mentioned  previously. 

The  following  are  observed  from  the  equations  In  figure  9:  (1)  #imT 

varies  linearly  with  kf,  quadratlcally  with  Sf^  and  Inversely  with  Efn;  (2) 
j£ailC  also  varies  linearly  with  kf,  quadratlcally  with  SpQ  (assuming  fiber 
compressive  fracture)  and  Inversely  with  Ef n ;  (3)  j£q,22T  decreases  nonllnearly 
with  kf.  Increases  quadratlcally  with  Smj,  decreases  Inversely  with  Em,  and  In¬ 
creases  nonllnearly  with  Increasing  ratio  (Em/Ef22);  (4)  #8.220  and  #H12S 
are  matrix  dominated;  and  (6)  the  matrix  dominated  Impact  energy  densities 
decrease  nonllnearly  with  Increasing  void  content.  Several  examples  below 
Illustrate  use  of  the  equations  In  figure  9. 

Example  6.1.  Calculate  the  longitudinal  tensile  Impact  energy  density 
for  an  AS/E  (IMHS)  unidirectional  composite  with  0.6  fiber  volume  ratio. 

Using  respective  properties  from  table  1  and  kf  =  0.6  In  equation  1, 
figure  9. 

jPtllT  =  0.6  x  350  0002/2  x  31  000  000  =  1185  (lb/sq  1n.)/cu  in. 

Example  6.2.  Calculate  the  longitudinal  compressive  Impact  energy  density 
for  an  AS/E  (IMHS)  unidirectional  composite  with  0.6  fiber  volume  ratio. 

Using  Sfc  0.9Sfj,  Sfy  =  350  ksl,  Efn  =  31  mpsl  and  kf  =  0.6  In 
equation  2,  figure  9. 

#ltnC  =  °-6  x  (°-9  x  350  000)2/2  x  31  000  000  =  960  (lb/sq  In./cu  In.) 

It  Is  Instructive  to  calculate  #^nc  assuming  delamlnatlon/shear  fracture 
mode  (example  3.2).  For  this  case  #jnc  =  370  (lb/sq  1n.)/cu  In.  or  a  decrease 
of  about  61  percent. 

Example  6.3.  Calculate  the  transverse  tensile  Impact  energy  density 
(#822T)  Eor  an  AS/E  (IMHS)  unidirectional  composite  with  0.6  fiber  volume 
ratio.  Using  respective  property  values  from  tables  1  and  2  and  kf  =  0.6 
In  equation  3,  figure  9  and  degrading  Smj  for  Incomplete  Interfaclal  bond 
(example  3.3) 

smT  =  {l.O  -  (4  (0.05)/(l  .0  -  0.6)  ir]1/2)  15  =  9  ksl 

#822T  ■  [1.0  -  (  vTTT  -  0.6)  (1.0  -  0. 5/2.0)  l2 

x  [1.0  -vCTO.O  -  0. 5/2.0)]  x  90002/2  x  500  000 

=  26  (lb/sq  1n.)/cu  In. 

This  value  Is  about  2  percent  of  the  longitudinal  tensile  (1185  (lb/sq  1n.)/cu 
In.,  example  6.1)  and  Illustrates  the  fragile  nature  of  unidirectional  com¬ 
posites  when  subjected  to  transverse  loads. 
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Example  6.4.  Calculate  the  transverse  compressive  Impact  energy  density 
for  an  AS/E  (IMHS)  unidirectional  composite  with  0.6  fiber  volume  ratio. 
Recall  that  Incomplete  bond  does  not  degrade  the  transverse  compressive 
behavior  (example  3.4).  Using  respective  property  values  from  tables  1 
and  2  and  kf  =  0.6  In  equation  (4),  figure  9 
#l22c  -  [1.0  -  (/O  -  0.6)  (1.0  -  0. 5/2.0) ]2 

x  [1.0-  0. 6(1-0. 5/2.0)]  x  35  0002/2  x  500  000 

=  388  (Ib/sq  1n.)/cu  In. 

which  Is  about  15  times  22T  Indicating  substantial  "tougher"  behavior 
In  transverse  compression.  Also  1s  about  the  same  as  that  for 

AllC  (370  (lb/sq/ln./cu.  In.)  calculated  by  assuming  delamlnatlon/shear 
compression  fracture  mode  (example  6.2).  This  Imples  that  longitudinal 
compression  and  transverse  compression  fractures  probably  occur  simul¬ 
taneously  during  normal  Impact. 

Example  6.5.  Calculate  the  Intralaminar  shear  energy  density  for  an  AS/E 
(IMHS)  unidirectional  composite  with  0.6  fiber  volume  ratio.  Using  re¬ 
spective  property  values  from  tables  1  and  2  and  kf  =  0.6  In  equation  (5), 
figure  9  and  degrading  Sms  for  Incomplete  Interfaclal  bond  (example  3.4) 
•%I2S  •  [1.0  -  (y^6  -  0.6)  (1.0  -  0. 185/2.0)  l2 

x  [1.0  -  v/O  (1  .0  -  0.185/2.0)  x  97002/2  x  185  000 

1  ,  =  53.5  (lb/sq  1n.)/cu  In. 

Sms  =  {1.0  -  [4  (0.02)/(1.0  -  0.6)  *]l/2}  x  13  =  9.7  ksl 

The  effects  of  voids  on  matrix  dominated  Impact  energy  densities  can  be 
calculated  by  degrading  the  matrix  strength  first  using  equation  6,  figure  9 
and  then  substituting  the  degraded  value  In  the  appropriate  equation.  See 
also  example  3.6. 


UNIAXIAL  IMPACT  RESISTANCE— THROUGH-THE-THICKNESS 

Through-the-thlckness  Impact  resistance  In  unidirectional  composites  re¬ 
sult  from  out-of-plane  normal  Impacts.  These  are  defined  herein  as  Impact 
energy  densities,  are  denoted  by  the  generic  symbol  and  are,  respectively: 
(1)  longitudinal  Interlaminar  shear  (  iral3S) ;UT2)  transverse  Interlaminar  shear 
(#423S);  (3)  longitudinal-  flexure  (  Xw\r)  and  transverse  flexurq  (j£h22f)- 
Each  of  these  Impact  energy  densities  Is  associated  with  a  dynamic  stress  cor¬ 
responding,  respectively  to:  Sll3s,  Sa23S,  SlllF  and  Sa22F-  Tfiere  1s  also  a 
through-the-thlckness  normal  Impact  energy  density.  However,  this  Impact 
energy  density  Is  the  same  as  the  In-plane  Impact  energy  dens1tyj&22T  or  #o??r 
described  In  section  6. 

The  composite  micromechanics  equations  for  through-the-thlckness  Impact 
energy  densities  are  summarized  In  figure  10  with  attendant  schematics.  The 
following  are  observed  from  the  equations  In  figure  10:  (1)  X^-\2S  1s  the 

same  &ti2Sl  It  decreases  nonllnearly  with  Increasing  kf,  Increases  non- 
llnearly  with  Increasing  ratio  (Gm/Gf-|2) ,  It  Increases  quadratlcally  with  SmS 
and  Increases  Inversely  as  decreases;  (2)  %%23S  bas  about  the  same  behavior 
as  ^Jll3s;  (3)  -£jfiiF  Increase  linearly  with  kf  and  quadratlcally  with  SfT,  and 
Increases  Inversely  as  Ef11  and  the  ratio  (SfT/Sfc)  decrease;  (4)  jTo22F  de¬ 
creases  nonllnearly  with  Increasing  kf  and  with  increasing  ratio  Em/Ef22,  In¬ 
creases  quadratlcally  with  Smj,  Increases  Inversely  as  the  square  of  the 
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(SmT/Smc)  ratio  and  Increases  as  Em  decreases.  Several  examples  below  Illus¬ 
trate  use  of  the  equations  In  figure  10. 

Example  7.1.  Calculate  the  longitudinal  Interlaminar  shear  Impact  energy 
density  (  J£g.i 35)  for  an  AS/E  (IMHS)  unidirectional  composite  with  0.6 
fiber  volume  ratio.  Since  #9^3$  is  the  same  as  %i2S*  (the  equation  for 
4>113S  Is  Identical  to  eq.  5,  fig.  9)  from  example  6.5. 

•#113S  =  -#11 2S  =  53,5  (1b/sq.1n./cu.1n.) 

Example  7.2.  Calculate  the  transverse  Interlaminar  shear  Impact  energy 
density  ( ^i23S)  f°r  an  (IMHS)  unidirectional  composite  with  0.6 
fiber  volume  ratio.  Using  respective  property  values  from  tables  1  and  2, 
and  kf  =  0.6  In  equation  for  %%23S>  figure  10,  and  degrading  Sms  for 
Incomplete  Interfaclal  bond. 

smS  =  in*0  -  [4  ( 0 . 02 ) /( 1 . 0  -  0.6)  irl1/2i  13  =  9.7  ksl 
J£o23S  •  [1.0  -  0.6  (1.0  -  0.185/1.0) ]2  97002/2 

x  185  000  [1  -  0.6  (1.0  -  0.185/1.0]  =  67.7  (lb/sq  1n.)/cu  In. 

It  Is  Interesting  to  note  that  for  this  example  -#st23S  is  ab°ut  22  percent 
greater  than  4^135.  This  Increase  Is  mainly  due  to  Gf23  which  Is  about  50 
percent  of  Gf -j 2 *  Based  on  the  relative  values  for  ^9,135  and  %i23S  1nter_ 
lamlnar  damage  will  occur  first  due  to  dynamic  09,23- 

Example  7.3.  Calculate  the  longitudinal  flexural  Impact  energy  density 
(  i?9,i  1  p)  for  an  AS/E  (IMHS)  unidirectional  composite  with  0.6  fiber 
volume  ratio.  Using  respective  property  values  from  tables  1  and  2  and 
kf  =  0.6  In  the  equation  for.  figure  10  (assuming  SfC  =  0.9  Sfj) : 

jPallF  *  4.5  x  0.6  x  350  0002/31  000  000  x  (1.0  +  1.0/0.9)2 

=  2394  (lb/sq  In.)  cu  In. 

An  alternate  estimate  Is  to  use  S^u  =  210  ksl  from  example  3.1, 

%11C  =  11®  ^si  from  example  3.2  and  E^n  =  18.8  mpsl  from  example  5.1  In  the 
following  equation 

■%1F  =  4-5  s2l  1  T/e8.1  1  0.0  +  SallT/S9,11iC-)2 
j£lllF  =  4.5  x  210  0002/l 8  800  000  x  (1.0  +  210/118)2 

*  1366  (lb/sq.  In.)  cu.  In. 


It  Is  worth  noting  that  the  first  estimate  corresponds  to  S^nF  =  298  ksl 
for  fiber  compression  fracture;  the  second  estimate  corresponds  to  S^nF  =  227 
ksl  for  delamlnatlon/shear  compression  fracture  (example  4.4).  Also  the  second 
estimate  Is  about  57  percent  smaller  than  the  first  Indicating  that 
delamlnatlon/shear  Is  a  much  more  severe  fracture  mode  under  Impact. 


Example  7.4.  Calculate  the  transverse  flexural  Impact  energy  density 
( ^122F)  f°r  an  (IMHS)  unidirectional  composite  with  0.6  fiber  volume 
ratio.  Using  respective  property  values  from  tables  1  and  2,  degraded  Smj 
for  Incomplete  Interfaclal  bond  (eq.„(6),  fig.  9,  with  kv  =  0.05)  and 
kf  =  0.6  In  the  equation  for  J>i22F»  figure  10. 


15  =  9.0  ksl 


smT  a  J1.0  -  [4  ( 0 . 05 ) /( 1 .0  - 
^422F  =  4-5  [1.0  -  0.6  (1.0 

Fl.O  -  (  VO  -  0.6)  (1.0  -  0. 5/2.0) 


0.6)  ,r]l/2i 

-  0. 5/2.0)] 


1  +  9/35 


(9000/500  000) 
146~( lb/sq  In. )/cu  In. 
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It  Is  worth  noting  that  this  value  corresponds  to  Sfc22F  =  18.6  ks 1  which 
Is  the  lower  estimate  In  example  4.5.  Also  ^22F  1s  aboiJt  10  percent  of 
the  lower  estimate  In  example  7.3 

The  effect  of  voids  on  any  of  the  through-the-thlckness  Impact  energy  den¬ 
sities  Is  determined  by  degrading  Smy  or  Sm$  first  using  equation  6,  figure  9 
and  then  substituting  this  degraded  Sm  value  In  the  applicable  equation, 
figure  10.  The  remaining  steps  are  Identical  to  Examples  7.1  to  7.4. 


ENVIRONMENTAL  EFFECTS 

Environmental  effects  refer  to  the  effects  caused  by  the  presence  of  mois¬ 
ture  and  temperature  In  composites.  The  combined  effects  are  usually  called 
hygrothermal  effects.  Hygrothermal  effects  Influence  all  resin-dominated  prop¬ 
erties;  uniaxial  strengths,  fracture  toughness  and  Impact  resistance.  Hygro¬ 
thermal  effects  are  estimated  using  an  empirical  expression  (ref.  9).  The 
empirical  expression  and  the  rotation  used  are  summarized  In  figure  11.  Its 
application  to  uniaxial  composite  strengths  and  strength  related  properties  Is 
Illustrated  using  the  following  examples. 

Example  8.1.  Calculate  the  hygrothermal  effects  on  the  ply  transverse 
strength  assuming  AS/E  (IMHS)  unidirectional  composite,  kf  =  0.6, 

T  =  270°  F,  T  =  70°  F  and  1  percent  moisture  by  weight.  Several  steps 
are  required  for  this  calculation: 

1.  Tgd  =>  420°  F  table  2 

2.  Tgw  [0.005  (l)2  -  0.1  (1)  +  1.0]  420  =  380°  F 

3.  Hygrothermal  degradation  ratio  (Phtm/po)  f°r  resin-dominated  proper¬ 
ties  Is 

pHTM/po  =  [( 380-270) /( 420-70)  ]^2  =  0.56 

This  means  that  all  resin  dependent  properties  Em,  Gm  and  Sm  must  be  re¬ 
duced  by  this  ratio  prior  to  their  use  In  the  applicable  equation. 

4.  The  reduced  matrix  properties  (table  2)  are 
SmT  =  0.56  x  15  =  8.4  ks  1 

Em  =  0.56  x  0.5  =  0.28  mpsl 

5.  Degrade  Smj  for  partial  Interfaclal  bond  assuming  5  percent  voids  by 
volume 

SmT  =  M  -  [4.,(0.05)/(l  -  0.6)  k]\/2\qa  =  5.0  ksl 

6.  Using  equation  3,  figure  5  with  the  respective  degraded  properties 
Sj,22T  =  Cl  -P  -  (-y/O  -  0.6)  (1.0  -  0.25/2.0)]  5.0  =  4.2  ksl 

which  Is  a  decrease  of  46  percent  compared  to  room  temperature  dry  7.8  ksl 
(example  3.3).  Obviously,  this  Is  severe  degradation  of  the  hygrothermal 
environment  assumed  In  the  example.  It  Is  Important  to  note  that  the 
ratio  of  environmentally  degraded  to  room-temperature  dry  (4. 2/7. 8)  =  0.53 
which  Is  very  close  to  0.56  predicted  In  step  3  above.  This  Indicates 
that  the  hygrothermal  degradation  ratio  can  be  applied  to  either  (1)  resin 
or  (2)  resin-dominated  composite  properties  equally  well  (ref.  9).  It  Is 
recommended  to  use  the  glass  transition  temperature  of  the  composite  for 
the  second  case.  The  glass  transition  temperature  of  the  composite  Is 
about  50°  F  greater  than  that  of  the  resin. 
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Example  8.2.  Calculate  the  transverse  flexural  strength  for  the  same  com¬ 
posite  and  environmental  conditions  as  In  example  8.1.  Again  several 
steps  are  required 

1.  From  example  8.1,  step  3 
pHTM/po  =0.56 

from  which  follows:  Smc  =  0.56  x  35  =  1.9.6  ksl  and 

Em  =  0.28  mpsl 

2.  The  transverse  compression  stress  Is 

S 2,22c  =  [1.0  -  (  t/O-  -  0.6)  (1.0  -  0.28/2.0)]  19.6  =  16.7  ksl 

3.  Using  S222C  =  ^  anc*  %22T  =  ^sl  (exafnPle  8.1,  step  6)  In 

the  lower  estimate  equation  (example  4.5). 

S222F  -  3  x  4. 2/(1. 0  +  4.2/16.7)  =  10.1  ksl 

which  Is  about  0.54  of  18.6  ksl,  the  room  temperature  dry  value  calculated 
In  example  4.5.  This  calculation  also  Illustrates  that  the  hygrothermal 
degradation  can  be  applied  to  a  resin-dominated  composite  property. 

Example  8.3.  Calculate  the  transverse  fracture  toughness  Ce/222T)  F°r 
the  composite  and  hygrothermal  environment  In  example  8.1.  The  procedure 
for  this  calculation  Is  the  same  as  that  In  example  5.2. 

1.  Using  the  (Phtm/po)  =  8*56  In  the  equations  for  E222  and  8212* 
figure  8. 

E222  =  0 . 28/[  1 .0  -^0T  (1.0  -  0.28/2.0)]  =  0.839  mpsl 
6212  =  0.1 03/ [  1 .0  -y/O  (1.0  -  0.103/2.0)]  =  0.39  mpsl 

2.  E211  a°d  »2i2  remain  practically  unchanged. 

^211  =  anc*  WH2  =  0-26  (example  5.1,  steps  2  and  4) 

3.  %22T  =  ^sl  (examPTe  8-1*  step  6) 

4.  Substituting  respective  values  from  steps  1,  2,  and  3  In  equation  2, 
figure  7 

jdl 22T  -  4. 2/1.0  +  (0.84/18.7)1/2  x  [2  (1.0  -  0.26) 

+  18.7/0 . 39]1 /2  =  1 #7  ksl 

which  Is  about  0.53  of  the  value  calculated  In  example  5.2.  Even  In  this 
complex  expression  the  environment  degrades  the  composite  resin  dominated 
property  In  about  the'  same  ratio  as  the  resin  property. 

Example  8.4.  Calculate  the  transverse  Impact  energy  density  (  #222T)  for 
the  composite  and  environmental  conditions  In  example  8.1.  The  procedure 
for  this  calculation  Is  the  same  as  that  In  example  6.3. 

1.  The  degraded  properties  needed  for  equation  3,  figure  9  are 

Em  =  0.28  mpsl  (example  8.1,  step  4  and  SmT  =  5.0  ksl  (example  8.1, 
step  5) . 

2.  Using  respective  values  In  equation  3,  figure  9 
■%22T  =  [1.0  -  (  \/0'  -  0.6)  (1.0  -  0.28/2.0)  ]2 

x  [1.0  -y/iTF  (1.0  -  0.28/2.0)]  x  50002/2  x  280  000 

•  =  10.8  (lb/sq.  In.)  cu.  In. 

which  Is  about  42  percent  of  26  (lb/sq  1n.)/cu  In.,  In  the  room  tempera-: 
ture  dry  value  In  example  6.3.  This  ratio  corresponds  to  a  decrease  of 
about  equal  to  the  hygrothermal  degradation  ratio  raised  to  the  3/2  power 
or  (0.56)3/2.  indicating  again  that  resin  dominated  composite  properties 
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degrade  In  the  same  ratio  as  the  resin  when  subjected  to  hygrothermal 
environments . 


GENERAL  DISCUSSION 

The  several  examples  presented  Illustrate  the  usefulness  and  advantage  of 
having  a  unified  set  of  micromechanics  equations  summarized  In  Figs.  5  to  11 
for  the  strength,  fracture  toughness  and  Impact  resistance  of  composites.  The 
examples  also  Illustrate  how  the  various  strengths  and  other  mechanical  proper¬ 
ties  are  Interrelated.  In  addition,  they  provide  detailed  and  quantitative 
Insight  Into  the  micromechanic  strength  behavior  of  composites.  Furthermore, 
the  various  equations  can  be  selectively  used  to  conduct  parametric  studies  as 
well  as  sensitivity  analyses  to  assess  acceptable  ranges  of  various  constituent 
material  and  environmental  factors. 

Limited  comparisons  were  provided  between  predicted  values  and  available 
measured  data  for  some  of  the  numerical  examples.  It  Is  Important  to  note  that 
the  primary  purpose  of  this  report  Is  to  describe  a  unified  set  of  simple, 
working  equations  and  Illustrate  Its  versatility  with  a  variety  of  numerical 
examples.  These  examples  demonstrate  computational  effectiveness  and  Illus¬ 
trate  Interrelationships  of  various  strengths  and  other  properties  at  the 
mlcromechanlstlc  level.  It  Is  highly  recommended  that  the  reader  use  this 
unified  set  of  micromechanics  equations  to  predict  various  properties  of  Inter¬ 
est  to  him  and  compare  them  with  measured  data  or  with  known  values.  This 
provides  a  direct  approach  to  assess  the  application  and  limitations  of  these 
equations  as  well  as  guidelines  on  how  to  modify  them. 

Another  Important  aspect  of  having  this  unified  set  of  micromechanics 
strength  equations  Is  that  they  can  be  used  to  plan  and  guide  experimental 
programs  for  maximum  benefit  with  minimum  testing.  These  micromechanics  equa¬ 
tions  can  be  advantageous  In  a  number  of  other  ways.  Many  of  these  other  ways 
become  "self  evident"  after  some  familiarity  has  been  obtained. 

The  two  tables  summarizing  constituent  material  properties  Illustrate  the 
amount  of  data  needed  for  effective  use  of  a  unified  set  of  micromechanics 
equations.  The  data  In  these  tables  were  compiled  from  many  sources  and  many 
values  are  estimates  which  were  Inferred  from  predicted  results  and  curve  fits. 
The  data  are  Included  for  three  main  reasons:  (1)  to  Illustrate  that  the 
micromechanics  equations  need  numerous  properties;  (2)  to  bring  attention  to 
the  fact  that  many  of  the.se  properties  have  not  been  measured  and,  hopefully, 
to  stimulate  enough  Interest  to  develop  experimental  methods  to  measure  them; 
and  (3)  to  provide  Indicative  ranges  of  properties  of  both  fibers  and  matrices. 
It  cannot  be  overemphasized  that  the  data  should  be  considered  dynamic  In  the 
sense  that  they  should  be  continuously  modified  If  better  values  are  known  or 
become  available. 

Lastly,  the  unified  set  of  micromechanics  equations  described  herein.  In 
conjunction  with  classical  laminate  theories  and  combined  stress  failure  cri¬ 
teria  can  be  used  to.  calculate  laminate  strength  based  on  first  ply  failure. 
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CONCLUSIONS 


A  unified  set  of  composite  micromechanics  equations  of  simple  form  Is 
summarized  and  described.  This  unified  set  Includes  composite  micromechanics 
equations  for  predicting  (1)  ply  In-plane  uniaxial  strength;  (2)  through- 
the-thlckness  strength  (Interlaminar  and  flexural);  (3)  In-plane  fracture 
toughness;  (4)  In-plane  Impact  resistance;  and  (5)  through-the-thlckness  (In¬ 
terlaminar  and  flexural)  Impact  resistance.  Equations  are  also  Included  for 
predicting  the  hygrothermal  effects  on  strength,  fracture  toughness  and  Impact 
resistance.  Several  numerical  examples  are  worked  out  to  Illustrate  the  ease 
of  use  of  the  various  composite  micromechanics  equations.  The  numerical 
examples  were  selected.  In  part,  to  demonstrate  the  Interrelationship  of  the 
various  constituent  properties  In  composite  strength  and  strength  related 
behavior,  and  also  to  provide  comparisons  with  available  experimental  data. 
This  unified  set  of  micromechanics  equations  makes  It  possible  and  cost- 
effective  to  assess  composite  strength,  "fracture  toughness".  Impact  resist¬ 
ance  and  attendant  environmental  effects  for  preliminary  designs  of  composite 
structures. 
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TABLE  I.  -  FIBER  PROPERTIES3 


Name 

Symbol 

Units 

Boron 

HMS 

AS 

T300 

KEV 

S-G 

E-G 

Number  of  fibers/end 
Fiber  diameter 

Density 

Longit.  modulus 

Transv.  modulus 

Long,  shear* modulus 
Transv.  shear  modulus 
Long.  Poisson's  ratio 
Transv.  Poisson's  ratio 
Heat  capacity 

Long,  heat  cond. 

Transv.  heat  cond. 

Long.  th.  exp.  coef. 
Transv.  th.  exp.  coef. 
Long,  tensile  strength 
Long,  compression  str. 
Shear  strength 

"f 

df 

Of 

Efil 

Ef22 

Gf  12 
Gf23 
vf  1 2 
vf23 
°f 

KfH 

Kf22 

®fll 

“f22 

Sft 

sfc 

5fs 

in. 

lb/in.3 

106  psi 

10*  psi 

10f  psi 

106  psi 

btu/1 b/*  F 
btu/hr/ft2/’F/in. 
btu/hr/ft2/*F/in. 
10-f  In. /In. /  F 

10-6  in./in./’F 

ksi 

ksi 

ksi 

1 

0.0056 

0.095 

58 

58 

24.2 

24.2 

0.20 

0.20 

0.31 

22 

22 

•  2.8 
2.8 
600 
700 
100 

10  000 
0.0003 
0.070 
55.0 
0.90 
1.1 
0.7 
0.20 
0.25 
0.20 
580 
58 
-0.55 
5.6 
250 
200 

10  000 
0.0003 
0.063 
31.0 
2.0 
2.0 
1.0 
0.20 
0.25 
0.20 
580 
58 
-0.55 
5.6 
350 
260 

3000 

0.0003 

0.064 

32.0 

2.0 

1.3 

0.7 

0.20 

0.25 

0.22 

580 

58 

0.55 

5.6 

350 

300 

580 

0.00046 

0.053 

22 

0.6 

0.42 

0.22 

0.35 

0.35 

0.25 

1.7 

1.7 

-2.2 

30 

400 

75 

204 

0.00036 

0.090 

12.4 

12.4 

5.17 

5.17 

0.20 

0.20 

0.17 

21 

21 

2.8 

2.8 

600 

204 

0.00036 

0.090 

10.6 

10.6 

4.37 

4.37 

0.22 

0.22 

0.17 

7.5 

7.5 

2.8 

2.8 

400 

Transverse,  shear,  and  compression  properties  are  estimates  inferred  from  corresponding  composite  properties. 


TABLE  II.  -  MATRIX  PROPERTIES 


Name 

Symbol 

Units 

LM 

IMLS 

IMHS 

HM 

Poly 

PMR 

imide 

Density 

0m 

lb/ in. 3 

106  psi 

0.042 

0.046 

0.044 

0.045 

0.044 

0.044 

Modulus 

Shear  modulus 

Em 

Gm 

0.32 

0.50 

0.50 

0.75 

0.50 

0.47 

106  psi 

Poisson's  ratio 

0.43 

0.41 

0.35 

0.35 

0.35 

0.36 

Heat  capacity 

vc 

Btu/lb/'F 

0.25 

0.25 

0.25 

0.25 

0.25 

0.25 

Heat  conductivity 

Thermal  exp.  coef. 

Km 

Btu/hr/ft2/’F/in. 
10"°  inJin./"F 

1.25 

57 

1.25 

57 

1.25 

36 

1.25 

40 

1.25 

20 

1.25 

28 

Diffusivity 

Dm 

10" I®  in2/ sec 

0.6 

0.6 

0.6 

0.6 

0.6 

0.6 

Moisture  exp.  coef. 

0m 

in. /in. / M 

0.33 

0.33 

0.33 

0.33 

0.33 

0.33 

Tensile  strength 

Smt 

ksi 

8 

7 

15 

20 

15 

8 

Compression  strength 

TTC 

ksi 

15 

21 

35 

50 

30 

16 

Shear  strength  . 

Sms 

ksi 

8 

7 

13 

15  ' 

13 

8 

Tensile  fracture  strain 

cmt 

in. /in.  (  *  ) 

8.1 

1.4 

2.0 

2.0 

2.0 

2.0 

Compr.  fracture  strain 

in. /in.  (  X  ) 

15 

4.2 

5.0 

5.0 

4.0 

3.5 

Shear  fracture  strain 

Gms 

in. /in.  (  X) 

10 

3.2 

3.5 

4.0 

3.5 

5.0 

Air  heat  conductivity 

Btu/hr/ft2/*F/in. 

0.225 

0.225 

0.225 

0.225 

0.225 

0.225 

Glass  trans.  temp,  (dry) 

tgd 

F 

350 

420 

420 

420 

700 

700 

Notes:  LM  -  low  modulus;  IMIS  -  intermediate  modulus  low  strength;  IMHS  -  intermediate  modulus 
high  strength;  HM  -  high  modulus. 


Thermal,  hygral,  compression  and  shear  properties  aree  estimates  only;  =  Em/2  (1  +  vm). 


TABLE  III.  -  TYPICAL  EXPERIMENTAL  VALUES  FOR  UNIDIRECTIONAL 
COMPOSITE  (PLY)  UNIAXIAL  STRENGTHS*3*  (KSI) 


Fiber 

Epoxy 

Fiber 

volume 

ratio 

Longitudinal 

•  Transverse 

I n-p 1 ane 
shear- 

Flexural 

Inter- 
1 aminar 
shear 

tension 

compres. 

tension 

compres. 

*  Long, 

Trans. 

Boron 

505 

0.50 

230 

360 

9.1 

35.0 

12.0 

350 

40 

16.0 

AS 

3501 

.60 

210 

170 

7.0 

36.0 

9.0 

230 

18 

14.0 

hWS 

934 

.60 

120 

-90 

6.7 

28.5 

6.5 

T50 

(b  I7 

10.5 

T-300 

5208 

.60 

210 

210 

6.5 

36.0 

9.0 

260 

18 

14.0 

KEVLAR-49 

.60 

200 

40 

4.0 

9.4 

8.7 

90 

6 

7.0 

GLASS 

S  (901-S) 

1002S 

.60 

220 

'  120 

6.7  • 

25.0 

12.0 

320 

21 

14.0 

E 

/  ^  \ 

1002 

.60 

160 

90 

4.0 

20.0 

12.0 

,  165 

20 

.14.0 

*0*Data  in  this  table  was  compiled  from  refs.  2,  3,  4,  5. 
^'Estimate. 
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Figure  1.  -  Concepts,  math-models  and  equations  used  to  predict  unidirectional  composite  (ply) 
uniaxial  strengths  from  constituent  material  properties,  geometric  configuration,  fabrication 
process  variables  and  environmental  conditions. 


SHEAR  • 

Figure  2.  -Typical  fiber  composite  geometry. 
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Figure  3.  -  Typical  stress-strain  behavior  of  unidirectional  fiber  composites. 
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(c)  Transverse  tension.  (d)  Transverse  compression.  (e)  Intralaminar  shear. 

Figure  4.  -  In-plane  fracture  modes  of  unidirectional  (ply)  fiber  composites. 
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1.  LONGITUDINAL  TENSION: 

2.  LONGITUDINAL  COMPRESSION: 

FIBER  COMPRESSION: 
DELAMINATION/SHEAR: 

MICROBUCKLING: 
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Figures.  -  Composite  micromechanics:  uniaxial  strengths  -  in-plane. 
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Figure  6.  -  Composite  micromechanics:  uniaxial  strengths  -  through-the-thickne*ss. 
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Figure  7.  -  Composite  micromechanics:  uniaxial  fracture  "toughness". 
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Figure  8.  -  Composite  micromechanics,  mechanical  properties. 
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1.  LONGITUDINAL  TENSION: 

2.  LONGITUDINAL  COMPRESSION: 

3.  TRANSVERSE  TENSION: 

4.  TRANSVERSE  COMPRESSION: 

5.  INTRALAMINAR  SHEAR: 

6.  FOR  VOIDS: 


•%lT“kf  SfT/2Efll 
•%lC“kf  SfC/2Efll 
•%2T  «[l-(Vk^-Kf)  (1  -  Em/Ef22T]  2 

x[l  -Vkj  (1  -  Em/Ef22*]^mT^^m 

•%2C  “[l-^rM'l-EjM2 

x[l-Virf«-Em/Ef22)]S^c/2Em 

■*I12S  “D  -  (VT(  -  kf)  (1  -  Gm/Cfi2)  ]  2 

x[l-^f«l-Gm/Gfl2l]S2mS/2Gm 


Sm»|l-[4kv/fl,l(f)5r]1,2JSIn 


VOID- 


O  O  Ol 

o  o  ; 
,o  o  qi 


Figure  9.  -  Composite  micromechanics:  uniaxial  impact  resistance  -  in-plane.  (Energy  absorbed  per 
unit  volume  M^) 
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Figure  10.  -  Composite  micromechanics:  uniaxial  impact  resistance  -through-the-thickness, 
absorbed  per  unit  volume. ) 
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Figure  11.  -  Governing  equations:  micromechanics  -  hygrothermal  effects. 


